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Abstract 



Physical wave functions for the nucleon and the A + isobar are pre- 
sented, which unify the best features of previous models. With these 
wave functions we can calculate elastic form factors and the decays of 
the charmonium levels 3 Si, 3 Pi, 3 i"2 into pp in agreement with the data. 
A striking scaling behavior between R = |GmI/^m an< ^ the coefficient 
-E>4 of the Appell polynomial decomposition of the nucleon distribution 
amplitude is found; the implications for elastic nucleon cross sections are 
discussed. 

I. INTRODUCTION 

An interesting testing ground for applications of perturbative QCD has emerged 
in the study of exclusive processes and elastic form factors of few-quark systems. Such 
systems can be described within a convolution formalismEl assuming factorization of 
highly off-shell or large transverse momentum regions of phase space from regions of 
low momenta necessary to form bound states. Recent progress! in Sudakov-suppression 
techniques provides support for the conjectured infrared protection of the perturbative 
picture. 

For modelling of the nucleon and its low resonances, elastic form factors play a 
key role because they provide an integrated view of the implications of QCD from low 
to high Q 2 . Thus they offer a powerful link between theoretical concepts and measure- 
ments and they can serve to test both the scaling properties as well as the detailed 
structure of the nucleon wave function. The theoretical tools for such a description 
are provided by the hard-scattering amplitude which describes the perturbative quark- 
gluon interaction in a particular process, and the probability amplitude for finding the 
three-quark valence state in the scattered nucleon or nucleon resonance: $N(^i,Q 2 )- 
A major theme of this talk will be to examine how QCD deals with the derivation of 
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such distribution amplitudes for the nucleon and the A + isobar, focusing our attention 
on recent developments. 



II. GENERAL FEATURES 

The momentum-scale dependence of ^^(xi,Q 2 ) is given by 

M*UQ 2 ) =$as(x i )f;5 n <l n (x i )(^|) 7n , (1) 



n=0 



in which {$ n }^° are orthonormalized eigenfunctions of the interaction kernel of the 
evolution equation^ expressed in a truncated basis of Appell polynomials of maximum 
degree M, and $ as (xi) = is the asymptotic form of the nucleon distribution 

amplitude. The corresponding eigenvalues 7 n turn out! to be the anomalous dimensions 
of multiplicatively renormalizable J1/2 baryonic operators of twist three. Because the 
7 n are positive fractional numbers increasing with n, higher terms in this expansion are 
gradually suppressed. A basis including a total of 54 eigenfunctions (M = 9) together 
with the associated normalization coefficients and anomalous dimensions is given in 
Ref. |. 

The derivation of the nucleon distribution amplitude from QCD is intimately 
connected with confinement and employs nonperturbative methods. Using the proper- 
ties of the Appell polynomials, the inverse of Eq. (1) determines the (nonperturbative) 
expansion coefficients B n : 

N f 1 

B n (fi 2 ) = J^l idx]<S> n ( Xi )^(xufi 2 ), (2) 
so that the "renormalization-group improved" coefficients B n (Q 2 ) are given by 

asiQ2) da , .1 . 2 (\n(Q 2 /A 2 QCD ) 



In terms of the moments of the nucleon distribution amplitude, 

*£° j V)= [ [dx&x^Kixi,^), (4) 
Jo 



Eq. (3) becomes 



5 n (/i 2 ) _ VK ^ „ ^(ioj) 



/AT 120 ^ 



where the projection coefficients are calculable to any order M. Specifically, those 
up to order M = 9 have been tabulated in Refs. |,|5[ 

To determine the moments, a short- distance operator product expansion is per- 
formed at some spacelike momentum /1 2 where quark-hadron duality is valid.B One 



considers matrix elements of appropriate three-quark operators which are related to mo- 
ments of the covariant distribution amplitudesQ V, A, and T: $N(^i) — V( x i) — A{x\), 
$ N (1, 3, 2) + $ N (2, 3, 1) = 2T(1, 2, 3) with V(1,2,3)=V(2,1,3), A(l,2,3)=-A(2,l,3), and 
T(1,2,3)=T(2,1,3). 

III. DISTRIBUTION AMPLITUDES OF THE NUCLEON AND THE A+ 

ISOBAR 

Based on QCD sum-rule calculations, useful theoretical constraints on the mo- 
ments of baryon distribution amplitudes have been obtained .H&llll Physical wave func- 
tions and observables for the nucleorJli'i and the A + isobarffl are then calculated 
using the full set of these constraints to determine the first few expansion coefficients 
B Q in a truncated basis of Appell polynomials. Depending on the value of Aqcd, these 
models predict approximately the right size and Q 2 -evolution of G^, while they give 
^ = I^mI/^m — 0-5. An alternative nucleon distribution amplitude was proposed" 
to give |Crjvf| <C G P M , inaccordance to phenomeno logical data analyses^ and the lat- 
est high-Q 2 SLAC data" at the expense that some of the amplitude moments cannot 
match the sum-rule requirements! in the allowed saturation range.ll! 

However, several crucial questions have to be resolved: For instance, does the 
optimum solution to the sum rules automatically yield best agreement with the data? 
Do solutions exist with characteristics distinctive from those of the COZ and the GS 
amplitudes? If so, what are the fundamental ordering parameters to classify these 
solutions? In recent works00 we have shown that it is indeed possible to amalgamate 
the best features of COZ-typei and GS-typelH nucleon distribution amplitudes into a 
hybrid-like amplitude, we termed the "heterotic" solution (see Fig. 1). 

In order to develop a credible nucleon distribution amplitude, we employ a x 2 
criterion which parametrizes the deviations from the sum-rule intervals according to the 
moment order. This "hierarchical" treatment of the sum rules takes into account the 
higher stability of the lower-level momentsll! and does not overestimate the significance 
of the still unverified constraints! for the third-order moments. [For more details, see 
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The underlying assumption is that contributions of higher-order terms are either 
negligible or of minor importance relative to those of second-order. Then the model 
space is also truncated at states with bilinear correlations of fractional momenta and 
the pattern of solutions found in this order should dominate the (orthonormalized) 
Appell polynomial series at every order of truncation. In this way the parameter space 
of the Appell decomposition coefficients can be systematically scanned seeking for local 
minima of x 2 . Using for the first and second order moments either the COZ or the KS 
sum-rule constraints in conjunction with those of COZ for the third-order moments, a 
simple scaling relation between the ratio R and the expansion coefficient emerges 
as one progresses through the generated solutions.0 

We have plotted in the (£>4, R) plane interpolating solutions to the COZ sum 
rules (+ labels) and such to a combined set of KS/COZ sum rules (o labels). As it 
turns out (Fig. 2), there is no significant difference between the two treatments and 



this insensitivity justifies the whole approach. The presented curves are fits to the 
local minima of the COZ sum rules (solid line) and the KS/COZ sum rules (dotted 
line). They constitute an orbit with respect to x 2 i beginning in the heterotic region 
(small R and large positive B±) and terminating past the COZ cluster (large R and 
large negative B±). 

The lower part of the orbit is associated with the heterotic solution which cor- 
responds to the smallest possible ratio still compatible with the sum-rule constraints. 
The upper region of the orbit controls COZ-type amplitudes and contains a cluster of 
solutions densely populating the orbit in the R- interval 0.455 -j- 0.495 (see the inset in 
Fig. 2). This cluster contains the amplitudes COZ opt , KS/COZ opt which are associ- 
ated with the absolute minima of x 2 an d play the role of strange attractors for all other 
solutions with similar features.^ GS-type amplitudes correspond to local minima of x 2 
at considerably lower levels of accuracy and thus they constitute in the (B 4 , R) plane 
an isolated region (an "island") that is separated from the characteristic orbit by a 
large x 2 barrier. The profiles of the distribution amplitudes across the orbit change in 
an orderly sequence of gradations with some mixture of COZ and GS characteristics 
until the COZ amplitude is transmuted into the heterotic solution. El 



TABLE I. Theoretical parameters defining the nucleon distribution amplitudes discussed 



in the text. 


The "hybridity" 


angle $ is 


discussed 


in [19]. 










Model 


Bi 


B 2 


B a 


B 4 


B 5 


tf[dcg] 


R 


x 2 


Symbol 


Het 


3.4437 


1.5710 


4.5937 


29.3125 


-0.1250 


-1.89 


.104 


33.48 


• 


Het' 


4.3025 


1.5920 


1.9675 


-19.6580 


3.3531 


24.44 


.448 


30.63 


• 


COZ°^ 


3.5268 


1.4000 


2.8736 


-4.5227 


0.8002 


9.13 


.465 


4.49 


■ 


COZ u p 


3.2185 


1.4562 


2.8300 


-17.3400 


0.4700 


5.83 


.4881 


21.29 


+ 


COZ 


3.6750 


1.4840 


2.8980 


-6.6150 


1.0260 


10.16 


.474 


24.64 


□ 


cz 


4.3050 


1.9250 


2.2470 


-3.4650 


0.0180 


13.40 


.487 


250.07 


♦ 


KS lovj 


3.5818 


1.4702 


4.8831 


31.9906 


0.4313 


-0.93 


.0675 


36.27 


o 


KS/COZ^ 


3.4242 


1.3644 


3.0844 


-3.2656 


1.2750 


9.47 


.453 


5.66 


o 


KS up 


3.5935 


1.4184 


2.7864 


-13.3802 


2.0594 


13.82 


.482 


40.38 


o 


KS 


3.2550 


1.2950 


3.9690 


0.9450 


1.0260 


2.47 


.412 


116.35 


o 


Qgopt 


3.9501 


1.5273 


-4.8174 


3.4435 


8.7534 


80.87 


.095 


54.95 


A 


G5 mm 


3.9258 


1.4598 


-4.6816 


1.1898 


8.0123 


80.19 


.035 


54.11 


T 


GS 


4.1045 


2.0605 


-4.7173 


5.0202 


9.3014 


78.87 


.097 


270.82 


A 



Let us now turn to models with functional representations which make use of 
higher Appell polynomials in connection with additional ad-hoc cutoff-parameters.0H 
The inset in Fig. 2 shows how such models^ (stars) andil (light upside-down triangles) 
group around the optimum amplitudes COZ opt and KS/COZ opt , thus establishing the 
scaling relation between R and B4 in a much more general context. This result suggests 
that the inclusion of higher-order Appell polynomials in the nucleon distribution am- 
plitude is a marginal effect, as conjectured above. Those model amplitudes!!! which 
appear as isolated points scattered towards the GS island are unacceptable on physical 
grounds, either because they exhibit unrealistic large oscillations in the longitudinal 
momentum fractionsill or because they yield a wrong evolution behavior for the nucleon 
form factors." 

One place to test these results is in the data for the elastic cross sections cr p 
and <r n . For small scattering angles, where the terms oc tan 2 (9/2) can be neglected, 
there are two main possibilities for the ratio cr n /&- If the Dirac form factor F± is zero 
or small compared to the Pauli form factor Ffj^ then a n should be due only to the 
higher-order term F£. At large Q 2 the ratio would become (Mn is the nucleon mass) 

cT ^ (if) q2 an d would decrease with increasing Q 2 due to the extra power of 
1/Q 2 of the Pauli form factor. Alternatively, if F± is comparable to F£, then o~ n would 
eventually be due to F± at large Q 2 . Then the ratio cr n /cr p . would be given by some 

constant determined by the nucleon wave functions =>■ 0$) ■ m these expressions, 

the wave-function characteristics are parametrized by the (dimensionful) coefficients Cj, 
which are functions of the expansion coefficients B n and the "proton decay constant" 
|/ N | = (5.0 ±0.3) x 10- 3 GeV 2 . 

The principal result from the above discussion is that in the intermediate Q 2 
domain, a Q /a p should be within the range 0.238 and 0.01. Comparing with available 
datai, we see that the measured (T n /a p enters the estimated range already at Q 2 ~ 
8GeV 2 /c 2 (see Fig. 3). 

In view of these results, it is worth remarking that the present accuracy of QCD 
sum rules seems to be sufficient to limit cr n /a p within the observed region. Fig. 3 shows 
that the available data in the range Q 2 ~ (8~10)GeV 2 jc 2 are well below the calculated 
upper bound and still decreasing. This indicates that distribution amplitudes which 
give I Cm I / ~ 0-5 ma y be i n contradiction to experiment because they yield a Dirac 
form factor F° which starts to overestimate the data already at Q 2 ~ 8GeV 2 /c 2 . On 
the contrary, models which give a small value of |GmI/^m can ex pl arn the data only 
under the assumption that in this Q 2 region the Pauli contribution is still dominant. 
Fig. 4 serves not only to amplify the preceeding discussion but also to advertise the 
consistency of the heterotic model with the form factordata. A similar good agreement 
with the data is found also for the axial form factors.t9>y 

We now turn our attention to the A + isobar. It was pointed out0 that model 
amplitudes for the nucleon are characterized by an anticorrelation pattern between 
and G* M . COZ-like models yield |G^|/G^ < 0.5 and (GmI/^m small > while GS-like 



models lead to the reverse situation. This pattern was derived, under the assumption 
that the A amplitude can be crudely modelled by the symmetric part of the nucleon 
distribution amplitude. Fig. 5 shows that the transition form factor calculated with 
the nucleon heterotic amplitude and more realistic A amplitudes, derived from QCD 
sum rules,00 is positive with a magnitude between those of previous models. In or- 
der to obtain an optimum distribution amplitude for the A, we try to comply with 
the constraints of the CP0 and FZOZ0 analyses simultaneously.0 This concept leads 
to a hybrid-like amplitude, denoted again "heterotic". This solution fulfills all FZOZ 
constraints and provides the best possible compliance with the CP constraints. In 
addition, it gives the best agreement with the data (see Fig. 5 and Ref. |^). In par- 
ticular, when including the effect of perturbative (i.e., logarithmic) Q 2 evolution of the 
expansion coefficients B n , the combined use of the heterotic amplitudes for the nucleon 
and the A + yields a form factor behavior which conforms with the observed decrease 
of available data within their quoted errors.0 

There is yet another type of solution for the A amplitude — compatible with the 
sum rules00— but in sizeable disagreement with the data. This solution (FZOZ opt ) 
is obtained by demanding that calculated with COZ opt is positive. Thus, as in the 
nucleon case, optimum agreement with the (existing) sum rules does not automatically 
entail best agreement with the data. 

Exclusive decays of charmonium levels to pp are very sensitive to the nucleon 
distribution amplitude. The branching ratio for the decay of the Xci state ( J CP = 1 ++ ) 
into pp is proportional to the decay amplitude Mi, which involves <3>at and /n- Inputing 
the heterotic amplitude, Mi is computed using an elaborated integration routine which 
accounts for contributions near singularities.^ Thereby we find M^ ct = 99849.6 and as 
a result BR( 3 Pi — > ppj 3 P\ — ► all) = 0.77 x 10 _2 %, which is in excellent agreement 
with the recent high-precision experimental valuecB (0.78 ± 0.10 ± 0.11) x 10~ 2 of the 
E760 Colaboration at FNAL. 

Analogously for the Xc2 state (J PC = 2++), we find M 2 hct = 515491.2. Setting^ 
a s (m c ) = 0.210 ± 0.028, we then obtain BR( 3 P 2 -> pp/ 3 P 2 -> all) = 0.89 x 10" 2 % in 
excellent agreement with the FNAL valuefl (0.91 ± 0.08 ± 0.14) x 10" 2 %. 

Similar considerations apply also to the charmonium decay of the level 3 S\ 
with J PC = 1 . The partial width of J /if} (or Xco) into pp is r( 3 Si — > pp) = 



( na )6 1280 \U\ 
\' lLx s) 243tt M 



fa 



M 2 



Mq, where determines the value of the 3 Si-state wave function at 

the origin. Its value can be extracted from the leptonic width r( 3 Si — > e + e~) = (5.36 ± 
0.29)/ceV@ via the Van Royen- Weisskopf formula. The result is {mj/^ = 3096.93Mel/) 
\f\p\ = 409Me\A The heterotic amplitude gives for this transition M = 13726.8. Then, 
usingthe previous parameters, it follows that r( 3 S'i — > pp) = 0.14keV. From experi- 
ment!! it is known that T(pp)/T tot = (2.16 ± 0.11) x 10~ 3 with r tot = (68 ± 10)fceV, 
so that r( 3 Si — > pp) = O.^/ceV^ in remarkable agreement with the model prediction. 
The corresponding branching ratio is BR( 3 Si — > pp/ 3 S\ — > all) = 1.62 x 10 -3 with 
Ttot = (85.5^ g) fee V. To effect the quality of these predictions, we quote the results for 



the COZ amplitude:^ [Note that these authors use the rather arbitrary value a s = O.3.] 
BR( 3 P 1 -> pp/ 3 P 1 -> all) = 0.50 x 10~ 2 %, PP( 3 P 2 -> pp/ 3 P 2 -> all) = 1.6 x 10" 2 %, 
and T( 3 5i -»■ pp) = 0.34A;eV. 

IV. SUMMARY AND CONCLUSIONS 

Given the apparent success of the heterotic in predicting a variety 

of observables such as magnetic and transition form factors and several branching ra- 
tios of exclusive decays of charmonium into pp, it is optimistic to believe that this 
approach — albeit approximative for a complete analytical understanding of the nu- 
cleon distribution amplitude — is sufficient of reproducing the observed phenomena. 
Since higher than order M = 3 expansion coefficients are unspecified by the present 
knowledge of QCD sum rules, the model does not depend on unconstrained (higher- 
order) parameters. While higher-order effects on the nucleon distribution amplitude 
itself are found to be large,0@ the agreement with the data is actually not improved.0 
This is also true for the optimized version of the COZ amplitudeii, we have derived, 
which represents the global minimum of % 2 . We emphasize that the (normalized) coef- 
ficients B n calculated via the central values of the 10 independent sum rules of Ref. |8| 
do not correspond to a solution with x 2 = 0. Although such a solution must exist, its 
determination is not a trivial task. Furthermore, at relatively large distances probed 
in present experiments, still uncalculable contributions of higher twists are presumably 
more significant than higher-order terms of the Appell polynomial series. 

This work was supported in part by the Deutsche Forschungsgemeinschaft and the 
COSY-Julich project. 
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